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We present the first measurement of the B^-Bl oscillation frequency Am s . We use 1 fb _1 of 
data from pp collisions at ^fs = 1.96 TeV collected with the CDF II detector at the Fermilab 
Tevatron. The sample contains signals of 3,600 fully reconstructed hadronic B s decays and 37,000 
partially reconstructed semileptonic B s decays. We measure the probability as a function of proper 
decay time that the B s decays with the same, or opposite, flavor as the flavor at production, 
and we find a signal consistent with B®-B s oscillations. The probability that random fluctuations 
could produce a comparable signal is 0.2%. Under the hypothesis that the signal is due to B®-B° s 
oscillations, we measure Am s — 17.31 to' it (stat.) ±0.07 (syst.)ps -1 and determine |I4<i/Vts| = 
0.208 «'Z(cxp.)^o 8 6 (theo.). 

PACS numbers: 12.15.Ff, 12.15.Hh, 13.20.He, 13.25. Hw, 14.40.Nd 



Neutral B mesons (bq, with q = d,s for B d , B s ) oscil- a B q meson produced at proper time t = to decay as a 
late from particle to antiparticle due to flavor-changing 
weak interactions. The probability density P + (P_) for 
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B° q (B®) at time t is given by 

P±(t) = Y e ~ Fqt t 1 ± cos(Am^)] , 

where Am, is the mass difference between the two mass 
eigenstates B° H and B° L Q, and T q is the decay width, 
which is assumed to be equal for the two mass eigen- 
states. The mass differences Anid and Am s can be 
used to determine the fundamental parameters |Vtd| and 
|Vf a |, respectively, of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix [2j, which relates the quark mass eigen- 
states to the flavor eigenstates. This determination, how- 
ever, has large theoretical uncertainties. A measurement 
of Am s combined with Am d = 0.505 ± 0.005 ps" 1 US 
would determine the ratio |Vtd/Vt s | with a significantly 
smaller theoretical uncertainty, contributing to a strin- 
gent test of the unitarity of the CKM matrix. Ear- 
lier attempts to measure Am s have yielded a lower 
limit: Am s > 14.5 ps -1 H at the 95% confidence 
level (C.L.). Recently the D0 Collaboration reported 
17 ps -1 < Am s < 21 ps" 1 at 90% C.L. |(| using a large 
sample of semileptonic B s Q decays. 

In this letter we report a measurement of Am s us- 
ing data from 1 fb _1 of pp collisions at y/s = 1.96 TeV 
collected by the CDF II detector at the Fermilab Teva- 
tron. We begin by reconstructing B s decays in hadronic 
(B s — > Dfn~, D+7r~7r + 7r~) and semileptonic (B s — > 
Da^£~Pe,£ = e oi /j) decay modes using charged parti- 
cles only |8|. Using the method of maximum likelihood, 
we extract the value of Am s from the probability density 
functions (PDFs) that describe the measured time devel- 
opment of B B mesons that decay with the same or oppo- 
site flavor as their flavor at production. The proper decay 
time for each B s is calculated from the measured distance 
between the production and decay points, the measured 
momentum, and the B s mass tub, = 5.3696 GeV/c 2 
The B s flavor (b or b) at decay is determined unambigu- 
ously by the charges of the decay products. 

To identify the flavor of the B s at production, we use 
characteristics of b quark production and fragmentation 
in pp collisions. At the Tevatron, the dominant b quark 
production mechanisms produce bb pairs. The b and b 
are expected to fragment independently into hadrons. In 
a simple model of fragmentation, a b quark becomes a 
B s meson when some of the energy of the b quark is used 
to produce an ss quark pair. The b and the s bind to form 
a B s . The remaining s quark may form a K~ . Similarly, 
a b that becomes a B® is accompanied by a K + . One of 
the two techniques used to identify the production flavor 
of the B s is based on the charge of these kaons (same-side 
tag). The second technique uses the charge of the lep- 
ton from semileptonic decays or a momentum-weighted 
charge of the decay products of the second b hadron pro- 
duced in the collision (opposite-side tag). 

The hadronic and semileptonic decay modes are com- 



plementary. Due to the large branching ratio, the 
semileptonic decays provide a tenfold advantage in signal 
rate at the cost of significantly worsened decay-time res- 
olution due to the unmeasured v momentum. Semilep- 
tonic decays dominate the sensitivity to oscillations at 
lower values of Am s . The fully reconstructed hadronic 
B s decays have superior decay time resolution, and our 
large sample of these decays is the unique feature that 
makes CDF sensitive to much larger values of Am s than 
other experiments. 

The CDF II detector 9] consists of a magnetic spec- 
trometer surrounded by electromagnetic and hadronic 
calorimeters and muon detectors [10j . The key fea- 
tures for this measurement include precision vertex de- 
termination provided by the seven-layer double-sided in- 
ner silicon strip detector RTlL lial supplemented with a 
single-sided layer of silicon mounted directly on the 
beampipe at an average radius of 1.5 cm. The 96-layer 
outer drift chamber 14] is used for both precision track- 
ing and dE/dx particle identification. Time-of- flight 
(TOF) counters [15j located just outside the drift cham- 
ber are used to identify low momentum charged kaons. 

Charm and bottom hadrons are selected using a three- 
level trigger system that exploits the kinematics of pro- 
duction and decay, and the long lifetimes of D and 
B mesons. A crucial component of the trigger system for 
this measurement is the Silicon Vertex Trigger Jl6| , which 

selects events that contain B° s — > Dfir~ and D+7r~7r + 7r~ 
decays. The trigger configuration used to collect the 
heavy flavor data sample is described in |l7| . 

To reconstruct B° s candidates, we first select Df candi- 
dates. We use D+ -> ()>tt + , K*{892)°K+, and 7r+7r-7r+, 
with — ► K + K~ and K*° — + K + ir~; we require that <p 
and K*° candidates be consistent with the known masses 
and widths 0] of these two resonances. These Z?+ candi- 
dates are combined with one or three additional charged 

particles to form D+£~, Dfir~ or _D+7r~7r + 7r~ candi- 

i — o 
dates. The Dj and other decay products of a B s candi- 
date are constrained to originate from a common vertex 
in three dimensions. For the K*(892)°K + final state, 
we remove candidates that are consistent with the de- 
cay D + — ► K~ir + ir + . We use a likelihood technique to 
identify muons 0] and electrons (lij. 

Backgrounds are suppressed by imposing a require- 
ment on the minimum transverse momentum px 20] of 

the B s and by requiring that the B s and decay ver- 
tices are displaced significantly from the pp collision po- 
sition. We find signals of 3,600 hadronic B s decays and 
37,000 semileptonic B s decays. 

For the hadronic decays, the invariant mass distribu- 
tion has a signal centered close to tob s = 5.3696 GeV/c 2 
with a width of 14 to 20MeV/c 2 , depending on the decay 
mode. Candidates with masses greater than 5.5GeV/c 2 
are used to construct PDFs for combinatorial back- 
ground. To remove contributions from B s — > D* + n~, 
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B s — > Dfp~ and semileptonic and other partially recon- 
structed decays, we require the mass of the decay can- 
didates to be greater than 5.3 GeV/c 2 . For semileptonic 
decays we take into account several background contri- 
butions, including B meson decays to two charm mesons 
and real D s mesons associated with a false lepton. 

The decay time in the B s rest frame is t = 
kILtttiBs/pt], where Lt is the displacement of the B s 
decay vertex with respect to the primary vertex projected 
onto the B s transverse momentum vector. The factor k 
corrects for missing momentum in the semileptonic de- 
cays (k = 1 for hadronic decays). To improve the decay- 
time resolution, we use event-by-event primary vertex po- 
sition measurements when computing the B s vertex dis- 
placement. The signal decay-time distribution is modeled 
with P(ti,a ti ) = e(U) jT s e- r ^'G{t' - U,a u )dt\ where 
U is the measured decay time of the ith candidate, T s is 
the B s decay width, Q(x — n,a) is a Gaussian distribu- 
tion of the random variable x with mean \i and width a, 
and Oti is the estimated candidate decay-time resolution. 
The decay-time efficiency function e(t) describes trigger 
and selection biases on the decay-time distribution and is 
determined from Monte Carlo simulation. For semilep- 
tonic decays, the K distribution is determined from Monte 
Carlo simulation and is convoluted with the signal decay- 
time distribution. The missing transverse momentum 
from unreconstructed particles in the semileptonic de- 
cays is an important contribution to the decay-time res- 
olution. To reduce this contribution and make optimal 
use of the semileptonic decays, we determine the k dis- 
tribution as a function of the invariant mass of the D s £ 
pair, rriD a i- The r.m.s. width of the k distribution is 3% 
(20%) for m D , t = 5.2 GeV/c 2 (3.0 GeV/c 2 ). 

We estimate the decay-time resolution cr ti for each can- 
didate using the measured track parameters and their 
estimated uncertainties. We calibrate this estimate us- 
ing a large sample of prompt D + mesons |'2lt . which we 
combine with one or three charged particles from the 
primary vertex to mimic signal topologies. For hadronic 
decays the average decay-time resolution is 87 fs, which 
corresponds to one fifth of an oscillation period at the 
lower limit on Am s (14.5 ps _1 ). For semileptonic decays, 
the decay-time resolution is worse due to decay topol- 
ogy and the missing momentum of unreconstructed de- 
cay products. For example, at t = 0, at = 100 fs (200 fs) 
for mn e t — 5.2 GeV/c 2 (3.0 GeV/c 2 ) and increases to 
a t = 115 fs (380 fs) at t = 1.5 ps. 

The flavor of the B s at production is determined us- 
ing both opposite-side and same-side flavor tagging tech- 
niques. The effectiveness Q = eT> 2 of these techniques 
is quantified with an efficiency e, the fraction of signal 
candidates with a flavor tag, and a dilution T> = 1 — 2w, 
where w is the probability that the tag is incorrect. 

Opposite-side tags infer the production flavor of the B s 
from the decay products of the b hadron produced from 
the other b quark in the event. We use lepton (e and 



//) charge and jet charge as tags, building on tec hniq ues 
developed for a CDF Run I measurement of Am^ |22| . If 
both lepton and jet-charge tags are present, we use the 
lepton tag, which has a higher average dilution. 

The dilution of opposite-side flavor tags is expected to 
be independent of the type of B meson that produces 
the hadronic or semileptonic decay. The dilution is mea- 
sured in data using large samples of B~ , which do not 
change flavor, and B , which can be used after account- 
ing for their well-known oscillation frequency. The com- 
bined opposite-side tag effectiveness is Q — 1.5 ± 0.1 %, 
where the uncertainty is dominated by the statistics of 
the control samples. 

Same-side flavor tags are based on the charges of 
associated particles produced in the fragmentation of the 
b quark that produces the reconstructed B s . In the sim- 
plest picture of fragmentation, a n + (ir~) accompanies 
the formation of a B~ (B + ), a ir~ (7r + ) acccompanies a 

5° (B°), and a K~ (K+) accompanies a B° s (B°). In 
Run I, CDF established this method of production flavor 
identification in measurements of Amj l24l and the CP 
symmetry violating parameter sin(2/3) [25j. In this anal- 
ysis, we use dE/dx 19] and time-of- flight information in 
a combined particle identification likelihood to identify 
the kaons associated with B s production. Tracks close in 
phase space to the B s candidate are considered as same- 
side kaon tag candidates, and the track with the largest 
kaon likelihood is selected as the tagging track. 

— o 

The performance of the same-side kaon tag for B s is 
expected to be different than for B~ and B . We pre- 
dict the dilution using simulated data samples generated 
with the pythia Monte Carlo [2(|. Control samples of 
B~ and B are used to validate the predictions of the 
simulation. The effectiveness of this flavor tag increases 
with the p T of the B° s ; we find Q = 3.5% (4.0%) in the 
hadronic (semileptonic) decay sample. The fractional un- 
certainty on Q is approximately 25%. This uncertainty 
is dominated by the differences between data and simu- 
lation for kaons found close in phase space to the B s ji]} 
and for the performance of the same-side kaon tag when 
applied to B~ . 

If both a same-side tag and an opposite-side tag are 
present, we combine the information from both tags as- 
suming they are uncorrelated. The addition of the same- 
side kaon tag increases the effective sample statistics by 
more than a factor of three. 

We use an unbinned maximum likelihood fit to search 
for B s oscillations. The likelihood combines mass, decay 
time, decay-time resolution and flavor tagging informa- 
tion for each candidate, and includes terms for signal and 
each type of background. The fit is done in three stages. 
First, a combined mass and decay-time fit is performed 
to separate signal from background and to fix mass and 
decay-time models. Combined fits for B s mass and de- 



6 



cay width in hadronic samples and for decay width in 
the semileptonic samples yield measurements consistent 
with established values Next, flavor asymmetries are 
measured for background components. The last step is 
a fit for — B s oscillations; the mass and decay-time 
models and background asymmetries are fixed from the 
previous two stages. 

The signal PDF has the general form: 

S±{t l ,at l ,V i ) — 

e(u) J Y e ~ Fst ' [1 ± AVl C0S ( Am ^')] G(ti - 1', <*u ) dt ', 

where T>i is the ith candidate dilution, and ti, at iy Q, and 
e(t) have been defined previously. Following the method 
described in [2^, we fit for the oscillation amplitude A 
while fixing Am s to a probe value. When all detector 
effects (T>i , a ti ) are calibrated, the oscillation amplitude 
is expected to be consistent with A = 1 when the probe 
value is the true oscillation frequency, and consistent with 
A = when the probe value is far from the true oscilla- 
tion frequency. Figure^ (upper) shows the fitted value of 
the amplitude as a function of the oscillation frequency. 
The sensitivity of the measurement is defined by the max- 
imum value of Am s where A = 1 is excluded at 95% 
C.L. if the measured value of A were zero. Our sensitiv- 
ity is 25.8 ps _1 and exceeds the combined sensitivity of 
all previous experiments At Am s = f7.3 ps _1 , the 
observed amplitude A = f.03± 0.28(stat.) is consistent 
with unity, indicating that the data are compatible with 
B° s -B s oscillations with that frequency, while the ampli- 
tude is inconsistent with zero: A/cr^ = 3.7, where a a 
is the uncertainty on A. The negative amplitudes mea- 
sured at frequencies slightly below and slightly above the 
peak frequency are expected and are due to the finite 
range in signal decay time that is imposed by the trigger 
and selection criteria. The systematic uncertainty on A 
is mainly due to uncertainties on tr ti and T>i. Since the 
effect of these uncertainties on A and a a are correlated, 
the ratio Aj a .4 has negligible systematic uncertainty. 

The significance of the potential signal is evaluated 
from A = \og[C A=Q / L A=1 (Am s )], which is the loga- 
rithm of the ratio of likelihoods for the hypothesis of 
oscillations (A — 1) at the probe value and the hypoth- 
esis that A = 0, which is equivalent to random produc- 
tion flavor tags. Figure ^ (lower) shows A as a function 
of Am s . A minimal value of A — —6.75 is observed 
at Am s = 17.3 ps^ 1 . The significance of the signal is 
quantified by the probability that randomly tagged data 
would produce a value of A lower than —6.75 at any value 
of Am s . We repeat the fit 50,000 times with random tag- 
ging decisions, and we find this probability is 0.2%. 

Under the hypothesis that the signal is due to B®- 
B s oscillations, we fix A = 1 and fit for the oscilla- 
tion frequency. We find Am s = 17.31 to.il (stat.) ± 
0.07 (syst.)ps -1 and the range 17.01 ps -1 < Am s < 



17.84PS- 1 (16.96PS- 1 < Am s < f7.9fps" 1 ) at 90% 
(95%) C.L. All systematic uncertainties affecting A are 
unimportant for Am s . The only non- negligible sys- 
tematic uncertainty on Am s is from the uncertainty 
on the absolute scale of the decay-time measurement. 
Contributions to this uncertainty include biases in the 
primary-vertex reconstruction due to the presence of 
the opposite-side b hadron, uncertainties in the silicon- 
detector alignment, and biases in track fitting. The 
measured B®-B s oscillation frequency is used to derive 
the ratio \Vtd/V ts \ = ^/ s;^- As inputs we use 
m B o/m B o = 0.98390 [2!j with negligible uncertainty, 
Am d = 0.505 ± 0.005 ps" 1 and £ = 1.21 toml US- 
We find \Vu/V ts \ = 0.208 ^(exp.) ™(theo.). 

In conclusion, we present the first measurement of 
Am s . The precision of this measurement is better than 
2%. The value of Am s is consistent with standard model 
expectations |3ll| and with previous bounds. Our mea- 
sured value of Am s allows us to determine |Vtd/V( S | with 
unprecedented precision and can be used to improve con- 
straints on the unitarity of the CKM matrix and on sce- 
narios involving new physics. 
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